Neutrino-nucleus coherent scattering as a probe of neutron density distributions 
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Neutrino-nucleus coherent elastic scattering provides a theoretically appealing way to measure 
the neutron part of nuclear form factors. Using an expansion of form factors into moments, we 
show that neutrinos from stopped pions can probe not only the second moment of the form factor 
(the neutron radius) but also the fourth moment. Using simple Monte Carlo techniques for argon, 
germanium, and xenon detectors of 3.5 tonnes, 1.5 tonnes, and 300 kg, respectively, we show that 
the neutron radii can be found with an uncertainty of a few percent when near a neutrino flux of 
3 x 10 r neutrinos/cm 2 /s. If the normalization of the neutrino flux is known independently, one can 
determine the moments accurately enough to discriminate among the predictions of various nuclear 
energy functionals. 

PACS numbers: 25.30Pt 



I. INTRODUCTION 

The size of a nucleus is one of its most fundamental 
properties. Although the distributions of protons in nu- 
clei are well known, the neutron distributions are com- 
paratively poorly constrained. A precise measurement of 
neutron radii could have important implications in both 
nuclear physics and astrophysics. 

In nuclear physics the most common framework for 
predicting neutron densities is energy density functional 
theory (DFT) — Skyrme or relativistic — with param- 
eters that are at least in part fitted to other nuclear ob- 
servables or pseudo-data such as nuclear-matter proper- 
ties, root- mean-square radii, atomic masses, etc. [l|. Re- 
cent efforts to determine new generations of energy den- 
sity functionals have pointed to the complex correlations 
among the values of physical quantities that they predict 
[__,__(]. Among observables used to optimize functionals, 
the neutron form factor is particularly important because 
it determines the neutron skin and radii, which in turn 
are strongly correlated in density functionals with the 
symmetry energy and incomprcssibility of nuclear matter 
[4j. Precise measurements of neutron radii could, there- 
fore, significantly improve the predictive power of energy 
functionals. 

In astrophysics, the size of the neutron skin may have 
important implications for neutron stars. While there 
are accurate measurements of pulsar orbital periods and 
masses, the radius of neutron stars, as well as their mo- 
ments of inertia or gravitational redshift remain poorly 
constrained by observation, and must be provided by the- 
ory Such global properties as masses, radii, and 
composition, are determined by the equation of state 
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(EOS) of neutron-rich nuclear matter, which involves 
such quantities as the density dependence of the nuclear 
symmetry energy. As just noted, the symmetry energy is 
itself strongly correlated with the neutron skin Q. Pre- 
cise measurements on the neutron skin, therefore, pro- 
vide information about the equation of state of neutron 
matter and thus the size of neutron stars [HI ] . 

Traditional measurements of neutron radii have in- 
volved hadronic scattering, with typical uncertainties of 
order 10%. Parity-violating electron scattering, used by 
the PREX experiment at Jefferson Laboratory to mea- 
sure the radius of lead, is cleaner. The parity-violating 
asymmetry, i.e. the fractional difference in cross sec- 
tion between positive- and negative-helicity electrons, is 
roughly proportional to the weak form factor, which is 
the Fourier transform of the weak charge density. If 
one can measure the asymmetry (at a single Q 2 ) to 3% 
then one can determine the root mean square neutron ra- 
dius to 1%. The uncertainty on the neutron radius from 
PREX is about 2.5% 03 but ±1% may be possible in 
the future. 

The use of neutrino-nucleus coherent scatterin g t o 
probe the weak form factor was first proposed in Ref. [12J ■ 
The authors considered a one tonne 40 Ar detector, with a 
nucleus described by a simplified form factor. Their anal- 
ysis suggested that ton-scale detectors could replicate the 
10% uncertainty of hadronic scattering methods when 
used in conjunction with a source of Michel-spectrum 
neutrinos. While large detectors are required since neu- 
trinos interact only weakly, the theoretical interpretation 
of the results is straightforward and model independent. 
Neutrino-nucleus coherent scattering has been proposed 
for a number of other purposes as well, for example to 
detect supernova neutrinos [lH , to measure the Weinberg 
angle jl4l Il5j ]. to look for a neutrino magnetic moment, 
and to search for sterile neutrinos [lj|. In all these cases, 
a weak nuclear form factor must be either measured or 
assumed before useful information can be extracted. 

Both neutron and proton distributions in the nucleus 
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affect neutrino-nucfeus coherent scattering, but the neu- 
tron distribution has much more ieverage. In this paper 
we suggest the use of a Taylor expansion to write the 
nuclear-neutron form factor in terms of moments of the 
neutron density distribution. Using this expansion and a 
simple Monte Carlo simulation, we show that neutrino- 
nucleus coherent scattering can probe not only neutron 
radii, but also the higher-order moments of neutron dis- 
tributions. We use the examples of argon [l4|, germa- 
nium [I?} and xenon targets to show the expected ranges 
of sensitivity. 

The paper is organized as follows. In section [III we 
introduce the model used to estimate neutrino-nucleus 
scattering count rates, including in our discussion the 
Taylor expansion of the neutron form factor and the cal- 
culation of the moments of the neutron distribution in 
nuclear DFT. In section UTTl we present and discuss the 
results of the Monte-Carlo simulations. 



II. COHERENT SCATTERING AND THE 
FORM FACTOR 

We present in this section the details of the model, in- 
cluding the kinematics of neutrino-nucleus coherent scat- 
tering, the dependence of the neutron form factor on 
the moments of the neutron distributions, and the DFT- 
based calculations of the moments. 



A. Kinematics 

To calculate the cross section for neutrino-nucleus co- 
herent elastic scattering, we sum the contributions of 
each nucleon to the amplitude, which we then square 
and sum over available phase space. The resulting cross 
section, for spherical nuclei (neglecting small corrections 
from various sources) is [l8| 
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where E is the energy of the incoming neutrino, T is the 
nuclear recoil energy, M is the mass of the nucleus, Gp 
is the Fermi constant, and Qw = N — (1 — 4 sin 2 9w)Z 
is the weak charge of the nucleus (sin 2 6\y m 0.231). 
The cross section also contains the form factor F 2 (Q 2 ), 
which is a function of the momentum transfer (Q = 
2E 2 TM/(E 2 - ET)). One of the neglected corrections 
is from higher multipoles in odd-A nuclei, which include 
the effects of deformation as well as nonzero spin. The 
higher multipoles add terms to Eq. (JTJ) only at order Q 4 , 
and even those changes are much smaller than 0(1/ A) 
for the nuclei considered here. 



The form factor corrects for scattering that is not com- 
pletely coherent at higher energies. It encodes informa- 
tion about the nuclear densities through a Fourier trans- 
form, which in spherical nuclei takes the approximate 
form 
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where p n>v (r) are the neutron and proton densities. We 
have neglected effects due to the finite size of the nucle- 
ons, which alter the relation between the point-neutron 
density and the cross section by modifying the form fac- 
tor at high Q. These effects could easily be included and 
would barely change the results of our sensitivity analysis 
below. 

The separation of the neutron and proton terms in Eq. 
<j2j) makes it possible to write the form factor as 

F{Q 2 ) = -!- [F n (Q 2 ) - (1 - 4sin 2 9 w )F p (Q 2 )] . (3) 
Qw 

Since the coefficient of the proton form factor is small, 
the scattering depends mainly on the neutrons, making 
neutrino-nucleus coherent scattering well suited to mea- 
suring the neutron distribution. 

There are two primary types of neutrino sources to con- 
sider: neutrinos generated from fission processes in nu- 
clear reactors, and neutrinos from the decay of stopped 
pions. Reactor neutrinos have lower energy, resulting 
in correspondingly low nuclear-recoil energies. Because 
background can obscure low-energy recoil, we consider 
neutrinos produced from the decay of stopped pions. 
Stopped pions are produced in large quantities at both 
spallation sources and accelerator sources. An example 
of a spallation source is the Spallation Neutron Source 
at Oak Ridge National Laboratory, which hits a mercury 
target with a beam of protons. Pions arc produced, with 
negative pions captured in the target and positive pions 
coming to rest and decaying. The pions decay through 
7r + — > -\- p + . The muon neutrinos are monoenergetic 
with an energy of 29.9 McV. The muons then come to 

e+ + v e + Va- The prob- 



(1) rest and further decay via p^ 



ability that neutrinos v e or antincutrino v ^ are emitted 
in the range (E, E + dE) read 
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where is the mass of the muon. The energy of the 
neutrinos range up to ~ 52 MeV, which results in typical 
nuclear recoil energies on the order of tens of keV to 
100 keV. The momentum transfer associated with these 
energies runs up to ~ 100 McV. 
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To calculate the number of scattering events as a func- 
tion of recoil energy, we fold the neutrino spectra with 
the cross section: 

w {T) = N tC j E ^ T m- { E,T) d E, (5) 

where N t is the number of targets in the detector, C 
is the flux of neutrinos of a given flavor arriving at the 
detector, the normalized energy spectra f(E) includes 
all three types of neutrino produced in pion decay, and 
E min (T) = i(T + \/T 2 + 2TM) is the minimum energy 
a neutrino must have to cause a nuclear recoil at energy 
T. The upper bound of is the maximum energy for 

a neutrino produced from muon decay at rest. 



B. Form-factor expansion 
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FIG. 1. (Color online.) The factor F(Q 2 ) for 40 Ar predicted 
by the Skyrme functional SkM* (solid black curve), and trun- 
cations of the expanded form factor at various orders of Q: Q° 
(dashed blue curve), Q 2 (dashed red curve), Q 4 (solid green 
curve). Terminating the expansion at Q 4 (with coefficient 
■^{Rn}) gives good agreement with the full form factor over 
the range of Q 2 relevant for the scattering of neutrinos from 
stopped pion beams. 
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Written this way, the form factor is a sum of the even 
moments of the neutron density distribution. These mo- 
ments are straightforward to calculate from the density, 
and represent physically relevant and measurable quan- 
tities. Since the neutrinos we consider have relatively 
low energy, we can truncate the expansion after just two 
terms for lighter nuclei such as argon and germanium, 
and three terms for heavier nuclei like xenon. As an il- 
lustration, we show in Fig. [T]the theoretical neutron form 
factor predicted by the Skyrme functional SkM* [l9| for 
40 Ar. Including moments up to (i? 4 ) is sufficient to re- 
produce the full form factor curve over the relevant range 
of Q values. In other words, we can fit experimental scat- 
tering data in 40 Ar with just two parameters, (i?^) and 
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Since the form factor is included in the calculation of 
the number of events, nuclei with different density distri- 
butions will produce different recoil-energy distributions. 
The recoil distributions therefore provide a good test for 
models that predict the density. We can increase the use- 
fulness of the recoil distribution by expanding the form 
factor in Q. The dominant neutron piece can be repre- 
sented as 

F n (Q 2 )~ J P n(r) (l-^ 2 + |f- 4 - 

«;v(i-^<i£> + S£<i£>- 



FIG. 2. (Color online.) Top: Event rates in 40 Ar as a func- 
tion of recoil energy, with two different RMS neutron radii. 
The red (solid) curve represents predictions of the Skyrme 
functional SkM*, while the blue (dotted) curve represents the 
same for an RMS radius made 10% larger, as described in 
the text. The flux at the detector is taken to be 3 x 10 7 
neutrinos/cm 2 /sec per flavor. Bottom: The difference be- 
tween the two curves on top. 



— -r 6 + • • • ) r 2 dr Fig.[2]shows the effects on event rates in 40 Ar of chang- 
/ ing a single important measure of the density distribu- 
te 6 , „ 6 ■ \r\ tion, the root-mean-square (RMS) neutron radius (R^) 5 . 
7! x nl F" We produced the fi gure as follows: First, we calculated 
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TABLE I. Isotopes and abundances used in the calculations 
for germanium and xenon (22j . 



Isotope 


Abundance 


Isotope 


Abundance 


' u Ge 


0.205 


128 Xe 


0.0191 


72 Ge 


0.274 


129 Xe 


0.264 


73 Ge 


0.078 


130 Xe 


0.041 


74 Ge 


0.365 


131 Xe 


0.212 


76 Ge 


0.078 


132 Xe 


0.269 






134 Xe 


0.104 






136 Xe 


0.089 



event rates as a function of recoil energy with the func- 
tional SkM*, via the expansion in Eq. ©. Next, we did 
the same with the RMS neutron radius 10% larger and 
no changes to the other terms in the expansion. The re- 
sulting curves for 40 Ar, with an assumed neutrino flux of 
3 x 10 7 neutrinos/cm 2 / sec per flavor, appear in the up- 
per panel of Fig. [2l while the difference between them 
is plotted in the bottom panel. A 10% difference in 
the RMS neutron radius results in a difference of ~780 
events/tonne/year, which is about 1.2% of the total event 
rate over the entire energy range. This difference is con- 
centrated at a nuclear recoil energy of 30 keV. 

We obtain similar results in Ge and Xe, though there 
we must use effective moments, which we define in the 
next section, to average over different isotopes. In Ge, a 
10% difference in the effective RMS radius yields an inte- 
grated difference of ^8100 events/tonne/year, or about 
6% of the total, concentrated around a nuclear recoil of 
15 keV. In Xe the same change in effective RMS radius 
results in a difference of ^20200 events/tonne/year, or 
about 8% of the total, and is concentrated at a nuclear 
recoil energy of 8 keV. 

As Q 2 goes to zero, the form factor approaches F(0) = 
1, so the low energy portions of the event-rate curves 
converge and all the difference curves go to zero. At high 
energies, there are very few events, so both event-rate 
curves approach zero and the differences also approach 
zero. The larger elements have smaller total recoil-energy 
ranges. In all three cases, however, the largest differ- 
ence in the event-rate curves occurs at a recoil energy 
of about 20% of the maximum. We note that the high- 



est and lowest energies will most likely be excluded from 
an experimental analysis because of background at the 
detector, and that the effect of changing (i? 2 ) 5 is most 
prominent in the energy range that is experimentally ac- 
cessible to c ryog enic detectors such as CLEAN [20[ and 
CLEAR [l4|, HJ, even though more events occur at inac- 
cessibly low recoil. 



C. Effective Moments 

Naturally occurring argon is made at 99.6% of 40 Ar: 
the expansion in Eq. ^ is therefore sufficient to accu- 
rately compute neutrino-nucleus scattering. Germanium 
and xenon, however, both have multiple isotopes that oc- 
cur naturally with relatively high abundance (see table 
IJ). To account for that fact, we define effective moments 
for these elements as follows: 

To calculate the event curves for germanium and 
xenon, it is necessary to sum over all isotopes. There are 
several isotope-specific quantities in the cross section, in- 
cluding the mass, neutron number, and the moments of 
the neutron and proton distributions. In addition, each 
isotope will have a different number of atoms in the de- 
tector. To calculate the rate for all isotopes, we therefore 
use 
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where the coefficient N t in Eq. ([5]) is now replaced by 
A^tMdetector and a summation over isotopes i of the cross- 
section with the weights Xi/Mi. Here, Xi is the natural 
abundance of isotope i, Mi is the mass of that isotope, 
Na is Avogadro's number, and Mdetector is the total mass 
of the element (including all its isotopes) in the detector. 

Since the form factor appears squared in the cross sec- 
tion, the sum in Eq. ([8]) will have neutron terms, proton 
terms, and terms that include both neutron and proton 
moments (see Eq. ([3])). The definitions for the proton ef- 
fective moments will follow the same pattern as those for 
the neutrons, so we will concentrate only on the neutrons 
at this point. 

Some algebra and the definitions of the cross section 
and form factor allow us to write the sum in Eq. ([5]) as 
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where (M) = J2 t X t Mi and Q 2 = 2E 2 T(M)/(E 2 - ET). 
In the above equation, we have kept all terms that have 
(i? 2 ), (i? 4 j), or both. (For xenon, we use all additional 
terms that include (i?^).) In this expression several ef- 
fective (isotope-weighted) moments occur, two of each 
order. The two effective second moments, after normal- 
ization, are 



(-R£)eff,l = 



(Rl) 
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and the two effective fourth moments are 

J2 t ^N 2 M 2 (Rj)„ 

E.^NfM 2 
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The differences between the values of the two moments 
of each order is small, as can be seen for the Skyrme 
functional SkM* in table [HI In fact, these differences are 
smaller than the typical numerical uncertainties in DFT 
calculations coming from the truncation of the basis or 
the size of the mesh. We therefore make the approx- 
imation that the two second moments are equal (call- 
ing them (i? 2 ) e g) and that the two fourth moments are 
equal (calling them (i? 4 ) e ff-) There are also terms in 
the cross section that involve (-R 2 ) 2 , from which we can 
define ((i? 2 ) 2 ) e ff. Although mathematically ((i? 2 ) 2 ) e ff is 
not equal to ((i? 2 ) e ff) 2 , numerically they are very similar, 
as shown in table [TTJ After equating them and making the 
similar approximations just described, we can rewrite Eq. 
[5] in terms of effective moments as 
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TABLE II. Numerical values for the different effective mo- 
ments of germanium and xenon as well as the percent differ- 
ence between definitions. The definitions for the effective mo- 
ments are given in equations 1 101 131 The values of ({-Rn) 2 )eff 4 



are compared to those of (R 



2x1/2 





Ge 


Xe 


{Rl)l^ (fm) 


4.0495 


4.8664 


{Rl)% (fm) 


4.0505 


4.8668 


% Difference 


0.02 


0.009 


(fm) 

{Ri)% (fm) 


4.3765 


5.2064 


4.3774 


5.2068 


% Difference 


0.02 


0.009 


i(Rl) 2 )\L A (fm) 
% Difference 


4.0509 
0.001 


4.8670 
0.01 



The approximations leading to Eq. (|T4|) from Eq. (J9J) 
cause an error of 0.01% over the entire event curve. 

We define effective moments for the protons in the 
same way as for neutrons. The terms in the cross section 
that involve both neutron moments and proton moments 
can be calculated from these effective moments. Since the 
proton moments are known quite accurately from elec- 
tron scattering, we are finally able to represent the recoil 
distribution in terms of just two unknown parameters, 
(B%) e « and (i? 4 ) c ff, or three for xenon. 



D. Density Functional Theory Calculations of 
Moments 



Our Monte-Carlo simulations require the knowledge of 
the radii and moments (R^ ) of the neutron and proton 
distributions in both even-even and odd-even isotopes. 



(14) 



r 



In this work, we compute these quantities in DFT with 
Skyrme functionals. We use nine common parameteri- 
zations of the Skyrme functional (SkM* Qjl, SkI3 [H, 
SLy4 S3, SLy5 ||, SkX ||,HFB9 SkO [27|, UN- 
EDF0 and UNEDF1 0). These functionals are char- 
acterized by relatively different nuclear matter, surface, 
and deformation properties, and therefore provide a good 
"statistical" sample of Skyrme functionals. 

We model pairing correlations with a density- 
dependent delta pairing force with mixed volume-surface 
characteristics and the Lipkin-Nogami prescription to ap- 
proximate particle number projection. For each element, 
we fit the strength of the pairing force to odd-even mass 
differences according to the general procedure outlined in 
0: in 40 Ar for argon, 72 Ge for germanium and 130 Xe for 
xenon. Only parameterizations UNEDF0 and UNEDF1 
are accompanied by specific prescriptions for the pair- 
ing channel. We compute the ground-state in odd-mass 
isotopes by performing systematic blocking calculations: 
for a given odd-mass isotope, we consider all blocking 
configurations within 2 McV of the ground state of the 
neighboring even-even nucleus, and take the ground state 
of the odd isotope to be the configuration yielding the 
lowest energy with the correct spin. 

We carry out all these calculations with the latest ver- 
sion of the DFT solver HFBTHO The code solves 



the Skyrme Hartree-Fock-Bogoliubov equations in a har- 
monic oscillator basis with axial symmetry and, there- 
fore, time-reversal invariance. We perform all calcula- 
tions in a full HO basis of 20 shells with a basis frequency 
parameter defined by hu> = 1.2 x 41/A 1 / 3 . In argon and 
germanium isotopes the basis is spherical, while in xenon 
isotopes, the ground states of which are weakly deformed, 
its deformation is (3 2 = 0.3. Such characteristics ensure 
excellent convergence of the results [1^] . We do the block- 
ing calculations in the equal filling approximation [30l |. 
which agrees with the full blocking prescription to within 
100 keV or less H3. 
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III. RESULTS AND DISCUSSION 



A. Monte-Carlo simulations 

We use a simple Monte Carlo simulation to give an 
idea of how accurately the nuclear moments can be de- 
termined. We assume that a detector, filled with ei- 
ther 40 Ar, natural germanium, or natural xenon, expe- 
riences a flux from the decay of pions at rest of 3 x 10 7 
neutrinos/cm 2 /sec in each flavor for one year. The neu- 
trino production rate at the Spallation Neutron Source at 
Oak Ridge National Laboratory [3, DAE5ALUS jHH 
and the European Spallation Source [34l ] ranges from 
1 x 10 15 neutrinos/sec to 3.5 x 10 15 neutrinos/sec of each 
flavor. A flux of 3 x 10 7 neutrinos/cm 2 /sec corresponds 
to detectors placed approximately 16 m from the source 
at SNS, 18 m from the source at DAEDALUS, and 30 m 
from the source at ESS. 

We perform a Monte Carlo simulation that includes 
statistical error and uncertainty on the beam normaliza- 
tion. Another significant source of systematic error is the 
uncertainty in the detection efficiency. We discuss this 
in Sec. Ill B, and for the purposes of the Monte Carlo 
assume 100% detection efficiency. We also assume that 
lcptons and photons produced by charge-current and in- 
clastic neutral-current scattering can be detected and the 
corresponding events efficiently rejected as background. 

We calculate a recoil curve assuming that the true nu- 
clear density distributions arc given by the Skyrme model 
SkM* . We place the events into bins based on energy, and 
then add Gaussian-distributed statistical noise to each 
bin. We then take the general form factor from Eq. [14] 
and use \ 2 minimization to find the optimal values of 
(i? 2 ), (i? 4 ), or the effective moments, and the beam nor- 
malization. In the case of xenon, we use the same proce- 
dure to find the optimal value of (R®). Typical values of 
X 2 range from 0.5 to 10. 

We do separate sets of runs, some assuming that the 
normalization of the flux is determined by other means, 
and some allowing for an uncertainty of ±10% in the nor- 
malization. Because background is anticipated to be sub- 
stantial at high and low energies, we exclude the highest 
and lowest bins from the \ 2 minimization. The energy 
ranges included are 5-120 keV for 40 Ar, 5-70 keV for Ge, 
and 5-40 keV for Xe. All energy bins are 10 keV wide 
except for the lowest, which is 5 keV wide. Finally, we 
assign confidence levels to closed areas on an (i? 2 ) vs. 
(i? 4 ) plot by running the Monte Carlo many times and 
dividing the number of times the minimum-x 2 result falls 
in that area by the total number of runs. 

With this setup, we vary the size of the detector un- 
til the (i? 2 ) 1/2 or (i? 2 )^ 2 inside the 91% confidence re- 
gion vary by only about ±5% from the best values. For 
that level of precision, 3.5 tonnes of argon are necessary, 
1.5 tonnes of germanium, and 300 kg of xenon. The re- 
quired detector mass decreases with atomic size because 
the event rate increases roughly as A'' 2 . 

The sizes of current and proposed cryogenic detectors 



can give an idea of the feasibility of this measurement. 
In the case of argon, the existing ICARUS T600 de- 
tector contains 500 tons of LAr [351 ] , suggesting that a 
detector big enough to for measure the form factor is 
feasible. For germanium, existing dark-matter detectors 
such as CDMS II (3(| consist of a few kg of germanium. 
The TEXONO-CDEX program is currently using a 1 kg 
high-purity germanium detector for neutrino physics and 
dark matter searches 
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The MAJORANA [3J and 
GERDA [39| double-beta decay experiments will soon 
deploy about 40 kg of germanium enriched in 76 Ge. One 
proposed experiment, GEODM [13, H^], would be made 
up of 300 ~5 kg Ge crystals, making a total mass of ~ 
1.5 tonnes. Existing xenon detectors, such as XENON100 
[40j and LUX [4l|, are made up of on the order of a few 
hundred kg of xenon, approximately the amount required 
for a form factor measurement. A proposed experiment, 



the LUX-ZEPLIN project, will use 1.5 tonnes of Xe [41 
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FIG. 3. Confidence regions in the (R 2 n ) 1/2 ~ (i? 4 ) 1/4 plane 
for an argon detector of mass 3.5 tonnes. The curves enclose 
confidence regions of 40%, 91%, and 97%. The colored verti- 
cal band shows the experimental result reported for the RMS 
radius, obtained from argon-carbon scattering, in Ref. p^ j. 
and the black crosses are the predictions of some commonly 
used Skyrme functionals, including the functional SkM* that 
we use to generate the "data." Numerical results are sum- 
marized in tables 1 1 1 1 1 and IIVI Top panel: the neutrino flux is 
allowed to vary by ±10 %. Bottom panel: the flux is assumed 
to be known exactly. 

The results of the analysis appear in Figs. [3j HI and 
[SJ The closed curves correspond to 40% confidence, 91% 
confidence, and 97% confidence. As mentioned above, 
we considered two cases: one in which the normaliza- 
tion of the flux is allowed to vary (by ±10%), and a 
second in which the normalization is kept constant. The 
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FIG. 4. Same as Fig. [3] but for effective moments in germa- 
nium, and without experimental result. 



top panels of each figure show the results with the flux 
unconstrained within that 10% range, and the bottom 
panels show the same results with the assumption that 
the flux is known perfectly. The colored vertical band in 
Fig. [3] shows a model-dependent experimental result for 
the RMS radius, obtained from argon-carbon scattering 
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in Ref. [42J. There is a clear discrepancy between that 
result and the predictions of the 9 Skyrme functionals 
selected for this study, labeled by small crosses in Fig. [3] 
(the outlyer in Figures corresponds to the SkX func- 
tional of [25| , which predicts systematically smaller radii 
than other functionals). Those functionals include SkM*, 
the one we use to generate the "data." This discrepancy 
is mentioned by Ozawa et ai, but no explanation is of- 
fered. While we marginalize over (i?^) g for xenon, the 
quantity is poorly constrained and not included in the 
plot in Fig. [5l Numerical results at the 91% confidence 
level for the mean, minimum, and maximum of the (ef- 
fective) RMS neutron radius and fourth moment, (and 
sixth moment in xenon) appear in Tabs. Hill and Hvl 



B. Discussion 



There are several noticeable differences among figures 
[U HI and O The first is the dependence on (i? 4 ) 4 ^ 4 . 
For 40 Ar, the range of plausible values is large compared 
to the range of (i?^) 1//2 . As the nuclei grow in size to 

germanium and then xenon, the range of (R^)^ gets 
smaller. In germanium and xenon, it is comparable to 
the range of the effective RMS neutron radius. 

We can explain this difference by isolating the effects 
of the (effective) fourth moment on the recoil distribu- 
tions. In 40 Ar, a 10% change in the fourth moment re- 
sults in approximately 0.2% more events, as compared 
to the 1.2% difference with a 10% change in the RMS 



4\l/4 
n/eff 



of 



FIG. 5. Same as Fig. 0] but in xenon. 



radius. In comparison, the same change in (R 
10% results in 1.3% and 3% more events for germanium 
and xenon, respectively, as compared to 6% and 8% from 
a 10% change in (R„)l(^ . The nuclear recoil energy at 
which this difference is concentrated decreases and the 
importance of the fourth moment relative to the RMS 
radius increases as the nuclear mass increases. 

Our ability to learn about the nuclear quantities 
(R„)l(^ and (i? 4 ) 4 ^ 4 obviously improves if we can get 
an independent handle on the beam normalization. This 
can be seen clearly by comparing the top and bottom 
panels in Figs. 3, 4 and 5. The range in (i? 2 ) 4 ^ 2 in all 
three elements shrinks to ±2% at the 91% level in the 
bottom panels, where the beam normalization is known 
exactly. Likewise, the range in (i? 4 ) 4 ^ 4 decreases for all 
three elements. The effect is most dramatic for °Ar, 
where the uncertainty decreases to ~ ±10%, but it exists 
in both germanium and xenon as well. 

As mentioned above, we consider the systematic er- 
ror of the uncertainty in detection efficiency. In order to 
study the effect, we remove the statistical error and ran- 
domly distribute systematic errors, proportional to the 
number of events in the bin, at the level of 10%, 1% and 
0.1% in each bin in an uncorrelated manner. At the 10% 
level, detectors of the size considered here lose the ability 
to make any useful measurement of the neutron radius. 
When wc lower the uncertainty in detection efficiency to 
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TABLE III. Numerical results at the 91% confidence level for the 3.5 tonne 40 Ar detector, the 1.5 tones Ge detector, and 
the 300 kg Xe detector with L„ allowed to vary by ±10%. The first column contains the element, the second the moment or 
effective moment considered in the corresponding row, the third the calculated values of the moments or effective moments for 
the Skyrme model SkM*, the fourth the mean values for the moments or effective moments, produced by the Monte Carlo, 
the fifth the percent difference between the mean values and the SkM* values, the sixth the minimum values chosen by the 
Monte Carlo, and the seventh the percent difference between the minimum and the mean value. The eighth column gives the 
maximum values chosen by the Monte Carlo, and the ninth column gives the percent difference between the maximum and the 
mean values. 







SkM* values 


Mean 


% Difference 
(from SkM*) 


Min 


% Difference 
(from mean) 


Max 


% Difference 
(from mean) 


40 Ar 


(Rl) 1/2 (fm) 
(Ri) 1 ' 4 (fm) 


3.4168 
3.7233 


3.4103 
3.6576 


-0.2 
-2 


3.2587 
2.8304 


-4 
-23 


3.5999 
4.3210 


+6 
+18 


Ge 


(Rl)l£ (fm) 
(Ri) 1 ^ 4 (fm) 


4.0495 
4.3765 


4.0516 
4.3603 


+0.05 
-0.4 


3.8792 
3.7276 


-4 
-15 


4.2697 
5.0096 


+5 
+15 


Xe 


(RDl's (fm) 
(RiVJs (fm) 
(Rl)l£ (fm) 


4.8664 
5.2064 
5.4887 


4.8648 
5.1914 
5.3149 


-0.001 
-0.3 
-3 


4.6788 
4.7180 
0.5491 


-4 
-10 
-90 


5.0980 
5.5521 
10.433 


+5 
+7 
+97 



TABLE IV. Same as table [TTT1 except for L„ held constant. 







SkM* values 


Mean 


% Difference 
(from SkM*) 


Min 


% Difference 
(from mean) 


Max 


% Difference 
(from mean) 


40 Ar 


(RlY r2 (fm) 
(Rt) 1/4 (fm) 


3.4168 
3.7233 


3.4154 
3.7018 


-0.04 
-0.6 


3.3483 
3.2826 


-2 
-11 


3.4933 
4.0865 


+2 
+10 


Ge 


(Rl) 1 ^ 2 (fm) 
(RifJs (fm) 


4.0495 
4.3765 


4.0491 
4.3679 


-0.009 
-0.2 


3.9857 
4.0826 


-1.6 
-7 


4.1175 
4.6546 


+1.7 

+7 


Xe 


(Rl)^ (fm) 
(Rit£ (fm) 
(Rl)l£ (fm) 


4.8664 
5.2064 
5.4887 


4.8654 
5.1990 
5.3877 


-0.02 
-0.14 
-1.8 


4.7958 
4.9265 
0.5491 


-1.4 
-5 
-90 


4.9323 
5.4478 
10.433 


+1.4 
+5 
+94 



1%, measurements of the neutron radius to ~ ±5-7% 
are possible. At this level, the range in the value of the 
fourth moment is ~ ±20%. If we lower the uncertainty 
in detector efficiency to 0.1%, the neutron radius can be 
measured to better than ~ ±1%, and the fourth moment 
to ±2%. 

IV. CONCLUSIONS 

Neutron radii are not only of fundamental interest for 
nuclear structure but are also needed to fully analyze 
supernova-neutrino signals (l3l | and interpret measure- 
ments of the Weinberg angle or of neutrino magnetic mo- 
ments (T3 |. At present the distributions of neutrons in 
nuclei is not known nearly as well as those of protons. 

We have presented a model-independent method, in- 
volving the Taylor expansion of the scattering form fac- 
tor, for extracting the RMS radius and fourth moment of 
the neutron density distribution in certain nuclei from the 
nuclear-recoil distribution in a neutrino-scattering exper- 
iment. The radius and fourth moment can also be be cal- 



culated theoretically, so that our technique will provide 
a straightforward connection between theory and exper- 
iment. To obtain a rough estimate of the effectiveness 
of this approach, we considered a stopped pion neutrino 
source of 3 x 10 7 neutrinos/cm 2 /s and liquid Ar, Ge, and 
Xe detectors in the tonne range. We conclude that it is 
possible to determine the neutron radius to a few percent 
if the uncorrclated error on the efficiency is less than 1%. 
The detailed analysis of the shape of the recoil spectrum 
in a cryogenic detector, such as the one we have suggested 
here, has not previously been considered. More detailed 
simulations of realistic experimental setups are therefore 
required for definitive feasibility studies. 
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